Coulomb excitation at intermediate energies 

H. Esbensen 

Physics Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 

(Dated: August 4, 2008) 

Abstract 

Straight line trajectories are commonly used in semi-classical calculations of the first-order 
Coulomb excitation cross section at intermediate energies, and simple corrections are often made 
for the distortion of the trajectories that is caused by the Coulomb field. These approximations 
are tested by comparing to numerical calculations that use exact Coulomb trajectories. In this 
paper a model is devised for including relativistic effects in the calculations. It converges at high 
energies towards the relativistic straight-line trajectory approximation and approaches the non- 
relativistic Coulomb trajectory calculation at low energies. The model is tested against a number 
of measurements and analyses that have been performed at beam energies between 30 and 70 
MeV/nucleon, primarily of quadrupole excitations. Remarkably good agreement is achieved with 
the previous analyses, and good agreement is also achieved in the few cases, where the B(EA) value 
is known from other methods. The magnitudes of the relativistic and Coulomb distortion effects 
are discussed. 

PACS numbers: PACS number(s): 23.20.-g; 25.70.De; 25.70.-z 
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I. INTRODUCTION 



The first-order excitation of a nucleus, induced by the Coulomb field from another nucleus, 
has been discussed in detail in Ref. The semi-classical perturbation theory developed 
there was based on a classical, non-relativistic Coulomb orbit for the relative motion of 
the projectile and target nuclei. A relativistic description was later developed but it was 
restricted to straight-line trajectories [2J. This leaves a gap in the theoretical description 
at intermediate energies (say, at 20-200 MeV/nucleon), where both relativistic effects and 
Coulomb distortions of the trajectory can be important. This is very unfortunate because 
many experiments have been performed in recent years at intermediate-energy, radioactive 
beam facilities. These include Coulomb dissociation experiments, for example, of 8 B, n Be 
and 11 Li, and Coulomb excitation measurements in inverse kinematics, mostly of the lowest 
2 + state in light and medium heavy nuclei j^]. 

The purpose of this paper is to devise a model that interpolates smoothly between the 
non-relativistic description of the Coulomb excitation which is based on Coulomb trajectories 
[l| and the relativistic description which is based on straight-line trajectories [2]. To develop 
an exact theory that contains the two extreme descriptions as limits is in general a difficult 
task. It has been studied by Bertulani et al. 

aa 

, who considered the effects of retardation 
explicitly for a Coulomb trajectory. The method proposed here is much simpler but it is 
sufficiently accurate for analyzing the data that have been obtained at intermediate energies. 
The point is that the experimental uncertainties are often large, typically of the order of 
10%, and there are also theoretical uncertainties, which can distort the analysis of data, for 
example, the influence of nuclear and higher-order processes. 

The accuracy of the proposed model of the intermediate-energy Coulomb excitation 
is tested by analyzing the results obtained with the commonly used 'Coulomb corrected' 
straight-line trajectory method, which originally was proposed by Winther and Alder 
Another test is to use B(E2) values that are known from other measurements (for example, 
of the life-time) and compare the calculated cross sections to measurements that have been 
performed at intermediate energies. 

There has recently been a debate in the literature about the validity of the analyses 
of Coulomb excitation experiments that have been performed in the past. It turned out 
to be caused by a misunderstanding of the experimental conditions, as pointed out in Ref. 
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However, independent of the controversy, it was claimed that corrections to the low- 
energy and high-energy theories of Coulomb excitation could be as large as 20% or 30% at 
intermediate energies [5|, [6|. In this paper it will be shown how this large uncertainty can 
be brought under control and reduced to only a few percent. 

The semi-classical theory of the non-relativistic Coulomb excitation is summarized in 
sections II to III. The expressions that are obtained from straight-line and Coulomb trajec- 
tories, respectively, are compared in section IV. The analytic expression for the relativistic 
Coulomb excitation in the straight line trajectory approximation is quoted in Section V, and 
it is used to devise a model which includes the combined effect of relativity and Coulomb 
trajectories. The model is tested in section VI against measurements and other calculations, 
and section VII contains the conclusions. 



II. NON-RELATIVISTIC COULOMB EXCITATION 

The semi-classical, non-relativistic description of Coulomb excitation [j| is summarized 
in the following. Thus we consider a target nucleus with atomic number Z 2 which is being 
excited by the Coulomb field from a projectile nucleus with atomic number Z\. It is assumed 
that the projectile and target do not overlap during the collision. We can therefore use the 
so-called far-field approximation which assumes that the intrinsic coordinate r of the target 
nucleus is smaller that the distance R(t) between projectile and target. The first-order 
amplitude for the electric excitation of the target nucleus, from an initial state \i > to a 

n 

final state \f >, is given by the multipole expansion 



where fwj = AEfi is the excitation energy, M AAt is the multipole operator r A V AA1 (f), and 



2A + IJ-00 Rit)^ 1 J-oo (R(t)) 

is the so-called orbital integral. It depends on the orbit R(t) of the projectile with respect 
to the target nucleus. The unit vector R = R/R is expressed in terms of the spherical 
coordinates (/?, 0) in the last expression. It is noted that this definition of the orbital 



integral, Eq. (j2j), differs by the factor w47r/(2A + 1) from the convention used in Ref. 
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The calculation of the orbital integrals is discussed in the following sections for two choices 
of the coordinate system. The first choice, system A, is the most convenient for numerical 
calculations that are based on a Coulomb orbit The second choice, referred to as system 
H, is more convenient at high energies, where a straight line trajectory becomes an accurate 
approximation and relativistic effects can be treated exactly . The transformation between 
the two representations will be discussed in order to be able to compare the two extreme 
methods and devise a model that interpolates smoothly between them. 



A. Cross sections 



The Coulomb excitation cross section is calculated as the product of the first-order exci- 
tation probability P/j and the elastic Rutherford cross section (da / dQ) r, 

The Rutherford cross section can be obtained from the Rutherford scattering formula, 
tan(6 l /2) = a/b, where b is the impact parameter and 

ZiZ 2 e 2 

a = iV' (4) 
is half the of distance of closest approach in a head-on collision, and M is the reduced mass, 
M = M 1 M 2 /(M 1 + Mi). The Rutherford cross section is 

/da\ 2nbdb 2 d / 1 \ 

WR = ^T = 7ra dn{ S m 2 (e/2))- (5) 

The average excitation probability obtained from Eq. (1) is the average over the initial 
magnetic sub-states and the sum over the final m-sub-states My, 

(2A + L)h p Hi + 1 M . Mf 

The last sum divided by (2ij + 1) is equal to (2A + times the multipole strength (or 
reduced transition probability) for the excitation, i. e., 
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III. ORBITAL INTEGRALS IN COORDINATE SYSTEM A 

In the coordinate system denoted by A (c.f. Ref. [lj]), the z-axis is perpendicular to the 
scattering plane so the orbit has the form R(t) = (x(t),y(t), 0). The angle (3 is fixed at 7r/2, 
and the orbital integral (2) can be written as 

S A - (-) r dt (x(t)+ ^ (t)) " (7) 

The factor d^ (ir/2) ensures that the orbital integral is non-zero only when A + /i is even. 

In the coordinate system A one chooses the x-axis as the symmetry axis of the Coulomb 
orbit so that x(—t) = x(t) and y{— t) = —y(t). From this it follows that = (S^) ^ s a 
real quantity. To calculate S£_„ when \i > one can use the expression 

S^-d^-)j_Jte {R[t)Y+ , +1 • (8) 

Below we discuss the calculation for a Coulomb trajectory and compare it to the result of 
the straight line approximation. 



A. Coulomb trajectories 

To evaluate the orbital integrals numerically for a Coulomb orbit one makes use of a 
dimensionless time- variable w (see Ref. so that 

R(t) = a[ecosh(w) + 1], t = ~[esinh(w) + w], (9) 

where a is defined in Eq. (j3J) and e is the eccentricity of the orbit, which can be expressed 
in terms of the impact parameter b or the scattering angle 9 in the center of mass system as 
follows, 



sin(0/3)' 

Inserting the Cartesian coordinates of the orbit (see 



x = a[cosh(u?) + e], y = aV e 2 — 1 sinh(w), z = 0, (11) 
into Eq. © one obtains, 

S i = ^<o(|) h„ (12) 
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where 

[°° [cosh(w) + e + iVe 2 - 1 sinh(w)]" 

J Am =/ rfw exp <« (esrnh w ) + w)\ , 13 

j-oo (ecosh(u>) + 1) A +^ 

and £ a = ua/v. Properties and even tabulations of I\n are given in Ref. jsj]. It will be 
calculated using a simple numerical integration with respect to w over the finite interval 
[-5:5], and using just a few thousand steps. The accuracy of the numerical integration can 
be tested in the case of a straight-line trajectory against the analytic expressions, which are 
discussed in appendix A. 



IV. COORDINATE SYSTEM H 



At high energies it is more convenient to use the coordinates system H where the z-axis 
is along the beam direction and the x-axis is along the impact parameter so that the y-axis 
is perpendicular to the scattering plane. The coordinates of the trajectory are therefore 
R(t) = (x(t),0,z(i)), which implies that Y Xfl (R) is real. Since Y£ = (— 1)^-^, it follows 
directly from the definition (2) that the orbital integral in coordinate system H must have 
the property, 

sf_» = (-irsf». (14) 

The coordinate system H is convenient at high energies because the analytic expressions 
for the orbital integrals that exist for straight-line trajectories are relatively simple in this 
representation even when the relativistic effects are included j^j. This feature will be ex- 
ploited in section V to devise an expression that contains the effects of relativity and the 
Coulomb distortion of the trajectory. To do that, we will need to transform the orbital 
integrals in coordinate system A to the new coordinate system H. 

To go from the H to the A coordinate system consists of a rotation of 7r/2 around the 
z-axis, followed by a rotation of —n/2 around the new y-axis, and finally a rotation of —n/2 
around the final z-axis. The required transformation is therefore 

sj£ = E<^<,(-f)s&- (is) 

n' 

Inserting Eq. (lT2"j) into this transformation we can write that 

= ^ where l£ = £ d^£) d£, (£) / V - (16) 
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Values of d*, (it/2) can be obtained from appendix D of Ref. [1]. The explicit expressions 



one obtains for dipole and quadrupole excitations are 
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(17) 
(18) 



A. Straight line trajectories 

In the straight-line trajectory approximation, the projectile moves with constant velocity 
v along the z-axis at an impact parameter b with respect to the target, R(t) = (6, 0, vt). The 
non-relativistic orbital integrals have the analytic form 

2 i x+l * 



X 



(A + / /)!(A-/i)! 



(19) 



where = ub/v is the adiabaticity parameter associated with the impact parameter b. The 
tilde on is a reminder of the straight-line trajectory approximation, and the superscript 
H refers to the coordinate system used here. It should be emphasized that the expression 
Eq. (fl9l) can be derived by inserting the straight-line approximation, Eq. (A3), into the 
transformation, Eq. (fT6|) . 



B. Coulomb trajectories 

In order to test the accuracy of Eq. (Tl9l) it is useful to write the general Coulomb 
trajectory result, Eq. ffl6l) . in a similar form, 



^ v /(A + / x)!(A-/i)! V ^ 



where 



Ktl'ib/a,^) = - - v/(A + fi)\(X - /*)! /J, 



1 / 1 \A 



(20) 



(21) 



and are the orbital integrals defined in coordinate system H according to Eq. ( TT6l) . They 
are given explicitly by Eqs. ( JT7j) and ( TTBl for dipole and quadrupole excitations. 

One can check that K e /J gives the correct modified Bessel function A" M (£) when one inserts 
the straight-line trajectory results in coordinate system A, Eqs. (A4-A6) of appendix A, 



into the definition f|T6|) of 1^. In the limit: £ a << 1 and b/a » 1, i. e., at high energies and 
large impact parameters, one should recover the result for a straight line, i. e. K^fJ (6/a, £ a ) 
— > K^b). This convergence will be illustrated in the next section. 
To evaluate the excitation probability we need the expression 

where 

« = ? (A + „)«*-/.)■ W ) ' <22) 

In the straight-line trajectory approximation, where K^fJ (b/ 'a, £ ) ~> K/j.(£b), one obtains 

^(6) = + and F 2 ^ b ) = ^ (3K^ b ) + 4K 2 ^ b ) + K 2 2 ^ b )) . (23) 

for dipole and quadrupole excitations, respectively. For — > these functions approach the 
constant values F\ — >1 and F 2 — >-l/3. 

C. Comparison of results 

Here the convergence of the Coulomb trajectories calculations towards the straight-line 
trajectory calculation is illustrated by comparing the functions F\(b/a,£ a ) defined in Eq. 
( l22l to the analytic expressions (l2"3i) for dipole and quadrupole excitations. The solid curves 
in Fig. 1 show the results of the straight-line trajectory approximation, Eq. f[2"3"j) . for dipole 
(A = 1) and quadrupole (A = 2) excitations, respectively, as functions of the adiabaticity 
parameter = cub/ v. The results of the Coulomb trajectory calculations defined in Eq. (|22|) 
are shown by dashed curves for three values of b/a. It is seen that the Coulomb trajectory 
results approach the straight-line result in a smooth manner for increasing values of b/a. 

The solid circles in Fig. 1 are the results one obtains by multiplying the straight-line 
trajectory expressions, Eq. (1231) . with the factor exp(— 7r£ a ) = exp(— ir£ b a/b) for b/a— 5. 
This correction factor was suggested by Winther and Alder Q and it is seen to repro- 
duce the dipole results (A=l) for the Coulomb trajectory with b/a=5 fairly well. How- 
ever, it does not work so well for quadrupole excitations (A=2). The problem is that the 
Coulomb trajectory results depend on b/a for — > 0, whereas the simple correction factor 
exp(— 7r£ a ) = exp(— irC, b a/b) is 1 in this limit. 
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The non-relativistic cross sections one obtains for the dipole excitation of the n Be l/2 + 
ground state to the l/2~ excited state and for the quadrupole excitation of 42 S to the 2 + 
state are illustrated in Fig. 2. The cross sections were calculated for a gold target with a 
minimum impact parameter of 14 fm and they are shown as functions of the beam energy. 
Although these cross sections are referred to as non-relativistic, it should be emphasized that 
the velocity that has been used here was actually obtained from the relativistic expression, 
Eq. (Bl), in terms of the beam energy. 

The excitation energies and B(EA) values that have been used for n Be and 42 S are shown 
in Table I, which will be discussed later. The top dashed curves in Fig. 2 show the straight- 
line trajectory calculations, and the thin dashed curves are the same results multiplied by the 
factor exp(— 7r^ a ). The solid curves are based on Coulomb trajectories. They approach the 
straight-line trajectory approximation fairly quickly for the dipole excitation but somewhat 
slower for the quadrupole excitation. It is also seen that the approximation of multiplying 
the straight line calculation with the factor exp(— 7r£ a ) works quite well for the dipole excita- 
tion when compared with the Coulomb trajectory calculation, whereas this approximation 
is poorer for quadrupole excitations. Other approximations have been applied to correct 
the straight line trajectory approximation for the distortion that is caused by a Coulomb 
trajectory and some of them will be discussed in section VI. D. 



V. RELATIVISTIC EXPRESSION 



The relativistic expressions for the orbital integrals in the straight-line trajectory approx- 
imation are Q], 



Sxp(rel) = - ' K^ b = -). (24) 

M 7^ v /(A + i u)!(A-/i)! KvJ l v 



where Gx^ can be extracted from Ref. 
the non-relativistic limit, whereas Ref. 
and quadrupole excitations one finds that 



2|. The notation used here is such that = 1 in 
2j uses a different normalization. For electric dipole 



G\q — G20 — G2±2 — —, G\±i — 1, C?2±i — -(1 H — n)- (25) 

7 2 7 Z 



It is seen that the relativistic effects on the Coulomb excitation enter into Eq. ( 1241) through 
the factor 7~ 1 Ga^ ) and in the adiabatic distance jv/cu of the adiabaticity parameter £ fe = 
u)b/(jv), which is the argument of the modified Bessel functions. 



To complete the discussion of relativistic effects one should also specify the kinematics 
of Coulomb scattering at relativistic energies. This is done in appendix B. One of the 
reasons large relativistic effects have been observed is actually due to the difference between 
relativistic and non-relativistic kinematics, the main one being the determination of the 
velocity from the beam energy, Eq. (Bl). A relativistic effect in Coulomb scattering is the 
7-factor which appears in the definition (B3) of half the distance of closest approach. This 

nn 

effect is commonly agreed upon [2|, [5( . There is also a relativistic correction to the reduced 
mass, Eq. (B2), which is less significant at intermediate energies, and it is usually ignored 



,2 



5]. However, for completeness, it is better to keep it in the case of really high energies. 
Finally, there are also relativistic effects in the transformation (B4) from the center-of-mass 
to laboratory scattering angles. 



A. Interpolating model 

To accurately calculate the Coulomb excitation at intermediate beam energies it is im- 
portant to include relativistic effects and the effect of the Coulomb distortion on the relative 
trajectory of projectile and target. It may be difficult to derive in a general expression for 
the Coulomb excitation amplitude because of the acceleration in a Coulomb orbit. However, 
one can devise a formula which gives the correct expression in the non-relativistic regime 
for a Coulomb orbit, and which reproduces the relativistic expressions for a straight line 
trajectory in the high-energy regime. We shall see that the following expression, 

serves the purpose of interpolating between the two energy regimes. 

It is first noted that Eq. (1261) is identical to Eq. (f20|) if we insert 7 = 1. The low-energy 

regime is therefore correctly described. In the high-energy regime, we can assume that 

b/a » 1 and £ a << 1, which implies that the straight-line trajectory limit will be reached, 

tt b 
K Xf/( b / a ^a) -> Kp(g b ), where & = 

according to the discussions in section IV.B and IV.C. Since the value of £ a in Eq. (l2"oT) is 
chosen as £ a = uia/^v) we obtain 

b uja ujb 

6 = - x — = — 

a 'yv 7f 
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which is the correct adiabaticity parameter for a straight-line trajectory in the relativistic 
regime, according to Eq. fl24|) . so the high-energy regime will also be described correctly. 



B. Total cross sections 

A great advantage of the straight-line trajectory approximation is that one can obtain 
analytic expressions for the total cross section, evaluated for all impact parameters larger 
that a certain minimum impact parameter bo . Thus one obtains [2J 

aA = *<WTWv) ( A) ? (a + fj.)\{\-fi)\ t) 9 ^ (27) 

where 

roo r -, 

g»(£) = 2vr j( Kl(£) = 7re 2 [^ + i(0^-i(0 - (28) 
according to Eq. (3.4) of Ref. [2j. For dipole and quadrupole excitations this yields 



<ta=i = 167r (f^) 2 [^2 So(£) + <?l(£)], (29) 

^= 2 = ^(f^) 2 i?(^2) (-) 2 [3p (0 +^i(07 2 (l + 7" 2 ) 2 + fla(e)]. (30) 

At intermediate and high energies, where one can assume that £& << 1, one obtains the 
following simple asymptotic expression for the quadrupole excitation cross section, 

This expression gives a surprisingly good estimate of the cross section at high energies and 
it provides a simple way of testing more elaborate numerical calculations. The expression 
(13~TT) shows that the high-energy cross section is insensitive to the excitation energy. One 
can also conclude that relativistic effects are not dramatic for E2 transitions because all of 
the 7 factors that appear in Eq. fl30l) disappear when one takes the limit — >0. 

We shall see in the next section that the large relativistic effects, which have been pointed 
out in the literature, are primarily caused by plotting the cross sections as a function of the 
beam energy T per nucleon. The cross section (13 ip . which is proportional to v~ 2 , will 



therefore be very sensitive to whether one uses non-relativistic (v = y2T/m) or relativistic 
kinematics (Eq. (Bl)) to determine the velocity. 
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For completeness it is noted that some relativistic effects do survive in the dipole cross 
section, Eq. fl29|) . when one takes the limit — >0 in the dipole cross section, Eq. fl29l) . 



This expression shows the well-know logarithmic dependence on 7. 

VI. APPLICATIONS 

Two examples of calculated cross sections are shown in Fig. 3, namely, for the dipole 
excitation of n Be to the low-lying l/2~ state, and the quadrupole excitation of 42 S to the 
lowest 2 + state, respectively. The cross sections are shown as functions of the beam energy. 
In both cases the minimum impact parameter was set to 60 = 14 fm and the target was Au. 
The excitation energy and multipole strength of the two transitions are given in Table I, 
which will be discussed below. 

The upper thick dashed curves in Fig. 3 show the relativistic Coulomb excitation cross 
section for straight-line trajectories. The solid curves are based on Coulomb trajectories 
and make use of the expression, Eq. fl26|) . for the interpolating model. The lower, thin 
dashed curves in Fig. 3 are the cross sections one obtains by inserting £ a = uoa/v into the 
expression, Eq. f l26l) . for the interpolating model. This is seen to be a poor approximation. 
Inserting instead £ a = ua/i^jv) (the thick solid curves) one obtains a smooth transition from 
the low-energy to the high-energy theory. 

A. Comparison to experiments 

A number of Coulomb excitation experiments have been performed at intermediate ener- 
gies and some of them are quoted in Table I. Those considered here are of interest because 
sufficient experimental information was provided to repeat the analysis. The measured cross 
sections, a exp , can be compared to the straight-line trajectory approximation, crstri, and to 
the cross section, o"cw, obtained from the interpolating model, Eq. It is seen that 

the latter model performs very well in comparison to the measured cross sections. The 
deviations are insignificant in comparison to the experimental uncertainty. 

It should be emphasized that most of the B(E2) values quoted in Table I were extracted 
from the data using the relativistic straight-line trajectory approximation with some correc- 




(32) 
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tions included for the Coulomb distortion of the trajectory. The good agreement between 
ocoui and the measured cross sections therefore shows that the previous analyses were very 
reasonable. The results of the straight-line trajectory approximation, astrii are quoted in 
the last column. They are typically 5 to 10% higher than the measured values. 

The ratio of measured and calculated cross sections is illustrated in Fig. 4. The solid 
circles show the ratio a exp /acoui with respect to the interpolating relativistic Coulomb ex- 
citation cross section. The average ratio is about 2% less than 1. The triangles in Fig. 4 
show the experimental ratio (T exp /astri to the relativistic straight-line trajectory calculation. 
Here the deviation from 1 is much larger. The deviation from the solid points reflects the 
significance of the Coulomb distortion of the trajectory. It amounts to about 3 - 9%. The 
largest deviations from the solid points occur in the low energy experiments, runs no. 4-8 
(see Table I.) 

There are two examples in Table I where the B(E2) values are known from other sources, 
namely, 24 Mg and 26 Mg. The calculated cross sections, acoui, are also in these cases consis- 
tent with the measured cross sections. This is very fortunate because the uncertainties are 
small (about 5% on the measured cross sections and 3% on the B(E2) values). These two 
measurements therefore provide an independent test of the interpolating model, Eq. (T261) . 

The last example on a comparison with data is the Coulomb excitation of 46 Ar [131 ]. 
which was measured for a range of maximum acceptance angles. The measurements are 
compared to two calculations in Fig. 5, namely, the straight-line and the 'interpolated' 
Coulomb trajectory calculation, Eq. ( 1261) . Both calculations are in good agreement with the 
data because the experimental uncertainties are much larger than the difference between the 
two calculations. 



B. Comparison to other methods 



The relativistic description for straight-line trajectories developed by Winther and Alder 
was based on the Lienard-Wiechert potential. An different method was used by Aleixo 
and Bertulani {4] who calculated the retardation effects explicitly for Coulomb trajectories. 
The latter approach was used by Bertulani et al. j^, (| to investigate the effects of relativity 
and Coulomb distortions at intermediate energies in much the same way it is done here. It 
is therefore of interest to compare Eq. (1261) to the predictions of the more elaborate theory. 
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An example of a comparison of calculated cross sections is shown in Fig. 6 for the Coulomb 
excitation of 46 Ar to the 2 + on a Pb target. The input is the same as used in producing 
table 2 of Ref. jf|, and it was assumed [l5] that the distance of closest approach for a 



Coulomb trajectory, r (b) = a + \/a 2 + b 2 according to Eqs. and (fTUj) . has the minimum 
value r m i n — R\ + R2 — 1.2(yl^ 3 + A^ 3 ). The solid line is the cross section obtained from 
the interpolating model, Eq. (|26|) . The upper dashed curve is the result for straight line 



trajectories with minimum impact parameter b min = \Jr m i n {r min — 2a). The lower dashed 
curve is the non-relativistic cross section for Coulomb trajectories. It is seen that the solid 
curve interpolates smoothly between the non-relativistic Coulomb trajectory calculation at 
low energy and the relativistic straight-line trajectory calculation at high energy. 

The solid points in Fig. 6A are based on the cross sections published in table 2 of Ref. . 
They have here been multiplied by the factor 1.44 because a factor of e 2 was unfortunately 
15 ] in the calculations of Ref. Q]. The results are presented in Fig. 6B in terms 



omitted 



of the ratio to the interpolating relativistic Coulomb trajectory calculation, Eq. ( l26l) . The 
average value of the solid points is close to one (actually 1.005 ± 0.002, to be precise) which 
shows that the two theories are in excellent agreement. This may not be so surprising because 
both theories approach the relativistic straight-line trajectory approximation at high energy 
and the non-relativistic Coulomb trajectory calculation at low energy. 

The lower dashed curve in Fig. 6B shows that relativistic effects are enormous at 500 
MeV/nucleon. It is interesting that the large relativistic effects have a very simple expla- 
nation. Thus, according to the asymptotic expression, Eq. ( 13T1) . the quadrupole excita- 
tion cross section is proportional to v ~ 2 at high energies. In the non-relativistic calcula- 
tion this implies cr^ R oc m/(2T). In the relativistic calculation one obtains instead of oc 
(m + T) 2 /[T(2m + T)], according to Eq. (Bl) of appendix B. The ratio of the two cross 
sections is therefore 

of* ^ m(2m + T) 

of ~ 2(m + T) 2 ' 1 ' 

(There is also a difference in the minimum impact parameter in the relativistic and non- 
relativistic calculations but the effect is small.) The expression ( 1331) is illustrated in Fig. 6B 
by the triangles which explain quite accurately the behavior of the non-relativistic calculation 
at high energy. 
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C. Significance of relativistic effects 

Another way to illustrate the effects of relativity is to recalculate the cross sections 
shown in Fig. 3 assuming that 7 = 1 everywhere in the underlying equations, Eqs. (I24H26I) . 
The velocity v, which appears explicitly in these equations, will be determined from the 
relativistic expression, Eq. (Bl), in order to avoid the trivial relativistic effect described by 
Eq. fl33|) . The results of such calculations show that the cross sections for the excitation 
of the 2 + state in 42 S, which were shown in Fig. 3, change by less than 0.5%, both for the 
straight-line trajectory calculation (124"]) and also for the interpolated Coulomb trajectory 
model, Eq. (j2Ej). 

Similar calculations performed for the dipole excitation of n Be (which were shown in Fig. 
3) change the cross section by less than 1% at energies below 200 MeV/nucleon. The change 
is about 10% at 1 GeV/u but that is not so surprising because the asymptotic dipole cross 
section, Eq. (132]) . does contain a 7 factor, whereas the asymptotic quadrupole cross section, 
Eq. (|3ll . does not. 

It is concluded that the relativistic effects in the Coulomb excitation of nuclei are large 
at intermediate and high energies but most of the effect is trivial and can easily be avoided 
by using the correct relativistic expression to determine the velocity from the beam energy. 
In the analysis of measurements it is also important to use relativistic kinematics when 
converting scattering angles into impact parameters. The non-trivial relativistic effects on 
the Coulomb quadrupole excitation, on the other hand, are surprisingly small. 



D. Significance of Coulomb distortion 

Let us finally take a look at how the straight-line trajectory approximation can be cor- 
rected for Coulomb distortions. Two examples are shown in Fig. 7, namely, the Coulomb 
excitation to the 2 + state of 16 N with excitation energy E x = 0.12 MeV (A), and also to the 
2 + state of 54 Ni with excitation energy E x = 1.4 MeV (B). The B(E2) values were taken 
from table 2 of Ref . |6( . The results are shown in terms of the ratio to the interpolating, rel- 
ativistic Coulomb excitation cross section, Eq. (|26|) . The calculations were again performed 

1 /3 1 /3 

with the minimum distance of closest approach r min = 1.2 (A-/ + A 2 ). 

The comparison with Ref. p] is illustrated in Fig. 7 by the solid points. The cross 



15 



sections from table 2 of Ref. [6| were again multiplied with the factor 1.44. It is seen that 
the solid points are very close to 1 for the excitation of 54 Ni. There are some discrepancies 
for 16 N, where the average ratio is 1.028 ± 0.012. It is therefore concluded that the two 
theories of relativistic Coulomb excitation agree within a few percent. 

The thick solid curves in Fig. 7 show the ratio for the relativistic straight-line trajec- 
tory approximation, and the dashed curves show the results of various ways to correct this 
approximation for Coulomb distortions. The top dashed curves show the straight-line ap- 
proximation multiplied by the factor exp(— 7r£ a ). This factor has a large effect for the heavier 
nucleus 54 Ni with the relatively large excitation energy but it has essentially no effect for 
the lighter nucleus 16 N with the small excitation energy. 

The lowest dashed curves in Fig. 7 show the straight-line trajectory result one obtains 
by replacing the minimum impact parameter b m i n by the effective value b e ff = b m i n + 
ixa/2. This approximation was justified in Ref. [2j for large impact parameters, where 
the excitation probability falls off exponentially as exp(— 2ujb/{^v)). Thus by multiplying 
the excitation probability for a straight line trajectory, Pstri, with the factor exp(— 7r£ a ) one 
obtains approximately, 

, . . , 2ub. . ua. r 2a;., ira 
Pstri exp(-7r£ a ) oc exp( ) exp(-7r — ) = exp (6+ — ) . 

ryv JV JV 2 

This argument does not always apply to the Coulomb excitation of low-lying states at 
intermediate energies because the minimum impact parameter is usually much smaller than 
the adiabatic distance '-/v/u. Using the effective minimum impact parameter has a very 
large effect on the calculated cross section. It produces a ratio in Fig. 7 that is almost as 
far below 1 as the ratio for the pure straight line trajectory calculation is above 1. The 
approximation is therefore not very useful for quadrupole excitations. It works apparently 
better for dipole excitations, as discussed in connection with Figs. 2 and 3, but that will 
not be discussed here. 

The second lowest curves in Fig. 7A and 7B are the results one obtains by choosing the 
effective minimum impact parameter, b e ff = r min , which is the minimum distance of closest 
approach that is used in the relativistic Coulomb trajectory calculation. This approximation 
was recommended by Goldberg [if]], and it is evidently the best choice out of the four 
examples of approximations shown in Fig. 7. 

In the analysis of an actual experiment the cross section ratio discussed above would 
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usually be closer to one because one would always choose a small acceptance angle (i. e., a 
large minimum impact parameter) in order to avoid the influence of nuclear and higher-order 
processes. 

VII. CONCLUSIONS 

A model has been devised for including relativistic effects in calculations of the Coulomb 
excitation cross section at intermediate energies. The model interpolates smoothly between 
the theory of non-relativistic Coulomb excitation at low energy and the relativistic, straight- 
line trajectory approximation high energy. The results that were obtained with this model 
compare very well with the calculations performed by Bertulani et al., who included the 
relativistic retardation effects explicitly for Coulomb trajectories. 

It was demonstrated that the large relativistic effects, which have been pointed out in 
the literature, are mainly caused by comparing calculations that are based on a relativistic 
and a non-relativistic velocity, respectively. The non-trivial relativistic effects, which are 
beyond the simple relativistic kinematics of two-body scattering, are surprisingly small for 
quadrupole excitations. 

The Coulomb distortion, which is responsible for the deviation between the straight line 
trajectory approximation and calculations that are based on a Coulomb trajectory, can 
also have a very large effect. However, the effects of the Coulomb distortion are usually 
suppressed by the experimental conditions and simple corrections can be made to improve 
the accuracy of the straight line trajectory approximation. 

The proposed interpolating model reproduces fairly well the analyses that have been 
performed previously, primarily of quadrupole excitation experiments, at beam energies in 
the range of 30 to 70 MeV/nucleon. The average deviation is only a few percent. The model 
also reproduces the measured cross sections in the few cases where the quadrupole excitation 
strength is known accurately from other sources. 

The good agreement with the previous analyses is partly due to the experimental con- 
ditions, which suppress the effects of the Coulomb distortion, and partly to the fact that 
some corrections for the distortion were made in the analyses. However, if high precision 
Coulomb excitation experiments were pursued, it would be necessary to treat the Coulomb 
distortion more accurately in the analysis. It is believed that the interpolating Coulomb 
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excitation model proposed here would provide a sufficiently accurate description. 
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VIII. APPENDIX A: STRAIGHT LINE TRAJECTORY 

In coordinate system A, a straight line trajectory has the coordinates x(t) = b, where 
b is the (constant) impact parameter, and y(t) = vt. Using the dimensionless integration 
variable s = vt/b, the orbital integral, Eq. (j7j), is therefore 

SY, = d x J-) — ds e^ bS 



AM fOl 2 ) vb x J ^ ^ + s 2)(A+ M +l)/2 

where = ub/v is the adiabaticity parameter. The tilde on is a reminder that we are 
using the straight-line approximation. An analytic expression for the integral is given in the 
book by Gragshteyn and Ryzhik Eqs. 8.432 no. 5, and one obtains 

^ = <4^j^Y + JkYeK.it>). w 



which is expressed in terms of modified Bessel functions of order n = (A + /x)/2. Here A + fi 
is even as mentioned earlier so n is an integer. We can express the result in a form similar 
to Eq. (JED with 

To evaluate this expression one can make use of the relations: 4-(x n K n (x)) = —x n K n _i(x). 
For dipole and quadrupole excitations one obtains 

Ji >±1 = 2& [k^tM^)}, (AA) 
ho = C [Jf 2 (&) " K (&] = 2 (^) 2 6, #i(&), (^5) 

4fc2 = ~^ 2 fe) T4^i(6) +3^0(6)], (^6) 

where the relation K%(x) = Kq(x) + 2fxK\(x) has been used in Eq. (A5). 

To improve the straight-line approximation, one can multiply the results, Eqs. (A4- 
A6), by the factor exp(— 7r£ a /2), according to Winther and Alder |2j. In fact, the exact 
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analytic expression, which has been obtained for a Coulomb trajectory in the case of dipole 
excitations, supports this suggestion; see Eq. (12) in appendix H of Ref. [l|. We will later 
on investigate how good this improvement and other approximations are for quadrupole 
excitations. 



IX. APPENDIX B: RELATIVISTIC KINEMATICS 

Here we specify the expressions that are used to calculate the relativistic Coulomb scat- 
tering. Most of them are taken from Jackson's book [181 ] . First of all, the kinetic energy T 
of the projectile is commonly given in units of MeV/nucleon so the 7 factor and the beam 
velocity v can be obtained from 



/T(2m + T) 

7m = m + T, j3 = v/c=- — , (Bl) 

m + T 

where m = 931.5 MeV is the nucleon mass (using the notation c=l.) 

The masses of projectile and target are denoted by Mi and M 2 , and the total energy in 
the center of mass system is (Jackson (12.31)) 



E' = J Ml + Mf + 2- MM. 



where 7 Mi is the laboratory energy of the projectile. Moreover, the energy and momentum 
of the projectile in the center of mass system are (Jackson, Eqs. (12.31-34)) 

Ml + £iM 2 Mi (Mi + 7 M 2 ) 



E< E' 
, M 2 MM 

where p = 7M1U is the momentum of the projectile in the laboratory frame. Note that the 
non-relativistic reduced mass MiM 2 /(Mi + M 2 ) has been replaced in the last expression by 
the relativistic reduced mass, 

MiM 2 M 1 M 2 , , 

M = — ^ = ■ (B2) 

E ^Ml + Mf + 2- M M, 
Rutherford's scattering formula in the center-of-mass system is derived from the trans- 
verse momentum transfer in elastic Coulomb scattering, estimated in the straight line ap- 
proximation by 

Ap± = p 1 sm(6) = — — — , 
19 



This is a reasonable approximation in high-energy forward-angle scattering but to match 
the non-relativistic expression, tan(#/2) = a/b, one should consider the Coulomb distortion 
of the trajectory. This would give a factor of cos 2 (#/2) on the right-hand side, so we obtain 

v ' ' p[vb -fM v 2 b 

Thus we recover the usual scattering formula, tan(#/2) = a/b, but the definition of a, Eq. 
01]), must be replaced by 

There are two corrections compared to Eq. (4). One is the factor I/7, which is commonly 
considered. The other is the relativistic reduced mass M , which is often replaced by the 
non-relativistic value. 

The scattering angle in the laboratory frame is determined by (see Jackson (12.50)) 

tan(0) = (54) 

7c.m.(cos(#) +a)' 

where 7 c . m . = (7M1 + M 2 )/E' (Jackson (12.35)), and 

= Mi Mi + 7 M 2 
Q M 2 1 M l + M2 

according to Jackson (12.54) for elastic scattering. It is seen that the transformation from 
the cm. to the laboratory system reduces to the usual expression for 7 — > 1. 
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TABLE I: Cross sections for the quadrupole excitation (A = 2) of different nuclei on a Au or Bi 
target, and the dipole excitation (A = 1) of n Be (last row). The experimental conditions are from 
the quoted references. The T is the beam energy per nucleon at mid-target and ^ (0S£L) is 
the maximum laboratory (center-of-mass) acceptance angle. The adopted B(E2) values [lo| for 
24 Mg and 26 Mg are also shown. The the last two columns show the calculated cross sections for 
the relativistic Coulomb (crcoui) and straight-line (crstri) trajectories. 



run 




E x 


B(EA) 


T 


xLab 
rvnax 




(JCoul 


O Strl 


no. 


Nucleus 


(MeV) 


(e 2 fm 2A ) 


(MeV / nucleon) 


(dee;) 


fmb) 


(mb) 


(mb) 


o 


24 Mg+Au [91 


1.3687 


467(28) 


36 


Qc.m. < 4 48 
"max — 


78.7(48) 


81.7 


88.0 




adopted [10] 




436(10) 








76.3(18) 


82.2 





26 Si+Au [9] 


1.7959 


336(33) 


41.8 


9%?- < 4.48 

max — 


55.8(55) 


56.3 


61.0 


1 


26 Me+Bi fill 


1.8087 


315 


67 


2.38 


44(2) 


45.9 


47.8 




adopted [101 




307(9) 








44.7(13) 


46.6 


2 


32 Mg+Au [11] 


0.885 


447 


71.2 


2.26 


91(10) 


92.7 


95.9 


3 


34 Mg+Bi [11] 


0.659 


541 


67 


2.38 


126(22) 


130.1 


134.9 


4 


38 S+Au [12] 


1.292 


235 


34.6 


4.10 


59(7) 


60.1 


64.7 


5 


40 S+Au [12] 


0.891 


334 


35.3 


4.10 


94(9) 


96.9 


103.4 


6 


42 S+Au [12] 


0.890 


397 


36.6 


4.10 


128(19) 


131.1 


139.9 


7 


44 Ar+Au [12] 


1.144 


345 


30.9 


4.10 


81(9) 


83.0 


89.6 


8 


46 Ar+Au [12] 


1.554 


196 


32.8 


4.10 


53(10) 


53.6 


58.3 


9 


46 Ar+Au [13] 


1.554 


212 


73.2 


2.90 


68(8 ) 


68.6 


71.9 




n Be+Au [14] 


0.32 


0.079 


57.6 


3.80 


244(31) 


246.1 


247.2 
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FIG. 1: The adiabaticity functions F\(b/a,£) are illustrated for A = 1 and 2. The top solid curves 
are the results for a straight line trajectory, Eq. (|23p . The dashed curves are obtained from 
Coulomb trajectories, Eq. (|22p . using the indicated values of b/a. The solid points were obtained 
by multiplying the straight line trajectory with exp{— 7r^ a ) = exp(— i^xib) f° r a/b=5. 



23 



1 



2000 



1000 



500 



1 



200 



2000 



1000 



500 



200 



100 



Straight line 
Strl*exp(-7t^ a ) 
Coulomb-traj 




Be on Au 
X=l 



10 



20 



Straight line 
Strl*exp(-7t^ a ) 
Coulomb-traj 



5 10 
T (MeV/u) 



20 



50 








50 



FIG. 2: Energy dependence of the non-relativistic cross sections for the dipole excitation of n Be 
and the quadrupole excitation of 42 S discussed in the text. The target is gold, and the minimum 
impact parameter was set to 14 fm in both cases. The upper thick dashed curves show the straight- 
line trajectory approximation, whereas the thin dashed curves have been corrected with the factor 
exp(— 7r£ a ). The solid curves are the results for Coulomb trajectories. 
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FIG. 3: Energy dependence of the relativistic cross sections for the dipole excitation of n Be 
and the quadrupole excitation of 42 S discussed in the text. The target is Au, and the minimum 
impact parameter was set to 14 fm in both cases. The thick dashed curves show the straight- 



line trajectory approximation, Eqs. (|29|30p . The solid curves show the interpolating, relativistic 
Coulomb excitation cross section, Eq. ([26]) . The thin dashed curves are also based on Eq. ([26]) 



but use £ a = uja/v in the second argument of Kt (b/a,£ a ). 
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FIG. 4: Ratio of measured and calculated Coulomb excitation cross sections for the different runs 
shown in Table I. Shown are the results for the interpolated relativistic Coulomb excitation (solid 
points) and the straight-line trajectory calculations (triangles). 
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FIG. 5: Cross section for the Coulomb excitation of the 2 + state in 46 Ar on a Au target at 73.2 
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FIG. 6: (A): cross sections for the 2 + excitation of 46 Ar on a Pb target. The minimum distance 
of closest approach was set to 11.41 fm. The relativistic straight-line (upper dashed curve) and 
Coulomb trajectory calculations (solid curve) are compared to the non-relativistic Coulomb trajec- 
tory calculation (lower dashed curve), and to the results of Ref. 6(] (solid points) which have been 
multiplied by 1.44. (B): ratios of the cross sections in (A) to the relativistic Coulomb trajectory 
calculation. The triangles show the estimate, Eq. f|33|) . 
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FIG. 7: Ratios of different cross sections to the relativistic Coulomb trajectory calculation. Results 
are shown for the excitation of the 2 + state in 16 N and 54 Ni, respectively, on a Pb target. The thick 
solid curves are the ratios for the straight line trajectory approximation. The dashed curves show 
various ways of correcting the straight-line approximation, such as multiplying with the factor 
exp(— TT^a), or using the effective minimum impact parameters discussed in the text. The solid 

n 

points are the ratios for the cross sections given in table 2 of Ref. [6j] multiplied by 1.44. 
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